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We investigate the metamagnetism in the periodic Anderson model with the fc-dependent c — / mixing 
by using the dynamical mean-field theory combined with the exact diagonalization method. It is found 
that both effects of the fc-dependent c — / mixing and strong correlation due to the Coulomb interaction 
between / electrons are significant for determining both the magnetization and the mass enhancement 
factor. For the case away from the half-filling, the results is consistent with the metamagnetic behavior 
observed in CeRu2Si2- 
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1. Introduction 

In the heavy fermion systems, the conduction (c) elec- 
trons mix with the almost localized / electrons, and form 
the strongly renormalized quasiparticles, which have effec- 
tive masses of 100^1000 times larger than the bare electron 
masses. Because of this strong renormalization, the heavy 
fermion quasiparticles often exhibit sensitive behavior to the 
external perturbation such as magnetic field H or pressure 
P. For example, the metamagnetism is one of the most re- 
markable features of the heavy fermion systems. The heavy 
fermion compound CeRu2Si2 exhibits an abrupt increase in 
the magnetization at an applied magnetic field H m ~ 7.7 T ' 
where the other physical quantities also show various anoma- 
lies without any symmetry breaking or magnetic phase transi- 
tion. This is known as a "heavy fermion metamagnetism" and 
has been studied extensively in past decades. '~ 5 ' 

To understand the physics of the heavy fermion systems 
including the metamagnetism from a microscopic viewpoint, 
we have to consider the effects of the c — / mixing and the 
strong correlation due to the Coulomb interaction between 
/ electrons which are crucial for the heavy fermion proper- 
ties. The Kondo lattice model (KLM) and periodic Ander- 
son model (PAM) are fundamental models for studying the 
heavy fermion properties and many theoretical studies have 
been carried out on the basis of these models. 6_9) Since it is 
difficult to solve these models exactly due to the many-body 
problem, some approximations have been required to solve 
these models, such as the Gutzwiller approximation, 6 ' slave 
boson mean-field, 7 ' non-crossing approximation 8 ' and 1/N 
expansion. 9 ' 

Dynamical mean-field theory (DMFT) has been known as 
one of the most powerful method to describe the strongly cor- 
related electron systems, since the DMFT includes the local 
quantum fluctuation effects which plays an essential role for 
the heavy fermion systems. 1(M5 ' The early DMFT studies on 
the PAM for the heavy fermion systems have succeeded in 
describing the strongly renormalized quasiparticles and the 
magnetization process. 12 ~ 14 ' These studies were restricted to 
the case with the half-filling and with a fc-independent c — / 
mixing, where the Kondo insulator state is realized. As for 



the metallic case, the PAM away from the half-filling has 
been studied by using the DMFT in combination with the it- 
erated perturbation theory (IPT). 15 ' In contrast to the exper- 
iment of CeRu2Si2, however, the metamagnetism occurs to- 
gether with the ferromagnetic transition in ref. 15). Hence, we 
need more realistic model for the study on the metamagnetism 
of CeRu 2 Si 2 . 

In addition to the strong correlation effect, the fc- 
dependence of c — / mixing is considered to be impor- 
tant for the metamagnetism of CeRu 2 Si2, 6,7 ' 9 ' where the fc- 
dependence originates from the local / state in the crystal 
electric field (CEF) 16 ' and yields the specific density of states 
(DOS) with the three peak structure. In this paper, we inves- 
tigate the heavy fermion metamagnetism on the basis of the 
PAM with the fc-dependent c — / mixing by using the DMFT 
combined with the ED (DMFT+ED) method. 

2. Model and formulation 

Our model Hamiltonian is the PAM with fc-dependent c — / 
mixing, which consists of the c electron term H c , f electron 
term Hf and the c — / mixing term H c f given by 



H — H c + Hf + H c f 
He = ^e k n c ka , 
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where ct is a creation operator for a c electron with wave 
vector fc and spin er =f, I and fj is that for a / elec- 
tron at site i with the Kramers doublet m = ± in the CEF 
groundstate. U is the Coulomb interaction between / elec- 
trons and Vkma is the c — / mixing. e k and ef m are energies 
for the c and the / electrons, respectively. For simplicity, the 
c electron is assumed to be a free electron with a rectangu- 
lar DOS and the band width is 2D. The CEF groundstate in 
the case of CeRu2Si2 is approximately given by J z — ±5/2 
Kramers doublet where the c — / mixing is given by V 2 (fc) 
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cos^fc = fc 2 /|fc|. 16) Here, we also assume that the magnetic 
field H is included only on the / electron e/ m = ej — mH 
with m = ± because the g value for / electrons is much 
larger than that for c electrons. 

We have solved the above Hamiltonian using the 
DMFT+ED method. In this method, the lattice Hamiltonian 
is mapped onto an effective impurity Anderson model, which 
includes the so-called Weiss field parameters (WFPs) and is 
solved with the Lanczos ED algorithm. The WFPs are deter- 
mined so as to satisfy the self-consistent equation given by, 

G m (ie v ) = — ( — e fm - £ TO (ie„) 
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where G m {ie v ) and H m (ie v ) are the Green's function and 
self-energy of the effective impurity Anderson model at 
T = 0, respectively. In eq. (5), ie v = i{2v + l)irf3 {y = 
0, 1, • ■ • ,n max ), and f3 is not an inverse temperature but 
a parameter which determines the energy resolution in the 
DMFT. n) 

In this study, we set the parameters as follows : the half of 
the c band width D = 2, c — / mixing strength V c f = 0.5, 
the frequency cutoff n max = 4000 and /3 = 600 - 800. We 
use the effective impurity Anderson model with the effective 
bath sites N s = 8 and confirm that the results are almost 
unchanged for N s = 6— 10. In the present paper, we will show 
the results for the half-filling case with (n) = (rif) + (n c ) = 
2.0 and the away from half-filling case with (n) = 1.7. 

3. Results 

Figure 1 shows the renormalized DOS p/(e) = 
— -ImG m (e + i0+) near the fermi level e = for several 
values of U at (n) = 2.0. In this calculation, we perform the 
analytic continuation of S m from the imaginary frequency to 
the real frequency by using the Padc approximation. When 
U = 0, we can see the three peak structure of the DOS which 
is due to the fc-dependent c — / mixing, where we refer the 
peaks in the left, center and right as peakl, peak2 and peak3 
as shown by arrows in Fig. 1, respectively. We note that there 
is no hybridization gap in the present model where the k- 
dependent c — f mixing vanishes around cos^fc ~ 1. When 
U increases, the renormalized quasiparticle band width de- 
creases with decreasing the integrated spectral weight of the 
quasiparticle band, and then the peakl and peak3 shift to the 
fermi level as shown in Fig. 1. In addition, we also observe 
broad peaks correspond to the Hubbard like bands around 
e/ = -C//2 ande/ + U = U/2 (not shown). 

To see the U -dependence of the correlation effect more ex- 
plicitly, we plot the renormalization factor defined by Z m = 
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(e) 
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for (n) = 2.0 and (n) = 1.7 in 

Fig. 2, where Z" 1 is nothing but the mass enhancement fac- 
tor: m*/m = Z^ 1 . In the cases with both of (n) = 2.0 and 
(n) = 1.7, Z m decreases with increasing U and the heavy 
fermion state is realized for U > 2.0, where Z" 1 becomes 
large as 10 or more. Here, we note that, for (n) — 2.0, we 
consider the particle-hole symmetry case with e/ = —U/2, 
while, for (n) = 1.7, e/ is fixed at tf = —1.7 measured rela- 
tive to the c band center. 

Fig. 3 (a) shows the magnetization curves M = (n^) — 
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Fig. 1 . The renormalized DOS near the fermi level e = for several values 
of U at (n) = 2.0 and e f = -U/2. 
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Fig. 2. The U dependence of the renormalization factor Z m for (n) = 2.0 
with t f = -U/2 (□) and (n) = 1.7 withe/ = -1.7 (A). 



(nj) for several values of U at (n) — 2.0. We can see nonlin- 
ear magnetization curves attributed to the fe-dependent c — / 
mixing. The metamagnetic behavior observed in CeRu 2 Si2 
is not obtained in the half-filling case. From the magneti- 
zation curve, we also calculate the differential susceptibility 
dM/dH (not shown) and find that dM /dH shows a peak at a 
critical magnetic field where the peakl and peak3 of the DOS 
cross the fermi level. When U increases, the critical magnetic 
field decreases with decreasing the quasiparticle band width 
shown in Fig.l. We note that the shape of the renormalized 
DOS are largely depends on H especially around the Hub- 
bard like bands which dominantly contribute to M for the 
case with large U where the / electrons are almost localized. 

The ^/-dependence of the mass enhancement factor Z" 1 
is plotted in Fig. 3 (b) for several values of U at (n) = 2.0, 
where Z + = Z- even for H ^ because of the the particle- 
hole symmetry. When H increases, Z" 1 decreases especially 
for the large U case, where the heavy fermion state is largely 
suppressed due to H. 

Fig. 4 (a) shows the magnetization curves for several values 
of U at (n) = 1.7. In this case, the fermi level with H = is 
located between the peakl and the peak2. When H increases, 
the differential susceptibility shows a peak at a critical mag- 
netic field H m where the peakl and the peak2 cross the fermi 
level almost at the same time. Then we observe the heavy 




fermion metamagnetic behavior at H = H m . In the interme- 
diate correlation regime with U ~ 1, the saturated magneti- 
zation above H m is almost a half of the full moment and such 
behavior is consistent with the result observed in CeRu2Si2- '' 
In the strong correlation regime with U ~ 4, the magnetiza- 
tion curve is similar to the case with the localized spin system 
where the metamagnetic behavior is suppressed. 

We also show the iJ-dependence of the mass enhancement 
factor Z- 1 at (n) = 1.7 for U = 0.8 and 2.4 in Fig. 4 (b). 
For both cases, Z+ 1 (majority component) is larger than ZZ 1 
(minority one) as previously obtained for the PAM with k- 
independent c — / mixing. 15 ' In addition, we also observe a 
shallow peak of ZZ 1 at H = H m due to the cooperative effect 
of the strong correlation and the fc-dependent c — / mixing. 

4. Summary and discussion 

In this paper, we have studied the heavy fermion metamag- 
netism on the basis of the PAM with the fc-dependent c — / 
mixing under the external magnetic field H, and have found 
that the both effects of the fc-dependent c — / mixing and the 
correlation due to the Coulomb interaction between / elec- 
trons are crucial for the metamagnetism. For (n) = 1.7, the 
results with the intermediate correlation is consistent with the 
metamagnetic behavior observed in CeRu2Si2- 

Experimentally, the heavy fermion metamagnetism is ac- 
companied by a large magnetostriction. 17 ' In the present 
study, we studied the case of the constant volume, i.e. the con- 
stant parameters efe, e/, V c f ■ When we vary the volume (pa- 
rameters) under a constant pressure, we may obtain such mag- 
netostriction together with an enhanced metamagnetic behav- 



ior. 9 ' The explicit calculation with the constant pressure to- 
gether with suitable parameter set of , e/ , V c / so as to repro- 
duce the experimental observation in CeRu2Si2 will be shown 
in a subsequent paper. 
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